function [Q]=fun_1(filepath2_DATA,filepath2_HEX,fileSave_Q,Area,step)

%% 读取仿真后数据
data1=load(filepath2_DATA);
data2=load(filepath2_HEX);
DATA = data1.DATA;
HEX = data2.HEX; 

%% 存储轨迹断裂点
[~,~,track_locat,~,~,~]=track(DATA,HEX);
%% 数据区域选择
if Area==1
    lat1=30;lat2=34;lon1=-99.5;lon2=-95.5;
end
if Area==2
    lat1=42;lat2=46;lon1=-116;lon2=-112;
end
if Area==3
    lat1=28;lat2=32;lon1=-101;lon2=-97;
end 
%% 计算每个区域的Q
[Q]=area_Q(DATA,HEX,track_locat,lat1,lat2,lon1,lon2,step); %步长a; %经纬度范围

save(fileSave_Q, 'Q'); %存储Q
end
